Mon. Not. R. Astron. Soc. 000,[T]{T5](2013) Printed 16 April 2013 (MN KTeX style file v2.2) 

Properties of the three-dimensional structure in the central region of 
the supernova remnant SNR 0540-69.3 * 



rn ; C. Sandin^t, P- Lundqvist^, N. Lundqvist^, C.-I. Bjomsson^, G. Olofsson^, 
^ Yu. A. Shibanov^'^ 

^ Leibniz-Institut fiir Astrophysik Potsdam (AIP), An der Sternwarte 16. 144 82 Potsdam, Germany 

^Department of Astronomy, Stockholm University, AlbaNova Science Center, Oskar Klein Centre, SE-106 91 Stockholm, Sweden 

^loffe Physical Technical Institute, Politekhnicheskaya 26, St. Petersburg, 194021, Russia 

"^St. Petersburg State Polytechnical Univ., Polytekhnicheskaya 29, St. Petersburg, 195251, Russia 



o 






o 

6 



> 

(N 

in 

o 

o 
m 



'i 



Submitted 12 Febi-uary 2013 



ABSTRACT 

We present and discuss new visual wavelength-range observations of the inner regions of the 
supernova remnant SNR 0540-69.3 that is located in the Large Magellanic Cloud (LMC). 
These observations provide us with more spatial and spectral information than were previ- 
ously available for this object. We use these data to create a detailed three-dimensional model 
of the remnant, assuming linear expansion of the ejecta. With the observations and the model 
we study the general three-dimensional structure of the remnant, and the influence of an ac- 
tive region in the remnant - a "blob" - that we address in previous papers. We used the fibre- 
fed integral-field spectrograph VIMOS at the Very Large Telescope of the European South- 
em Observatory. The observations provide us with three-dimensional data in [Oiii]/I5007 
and [Sii] /1/16717, 6731 at an 0'.'33x0'.'33 spatial sampling and a velocity resolution of about 
35kms ^ We decomposed the two, partially overlapping, sulphur lines and used them to 
calculate electron densities across the remnant at high signal-to-noise ratio. In our study we 
recover results of previous studies, but we are more importantly able to obtain more detailed 
information than before. Our analysis reveals a structure that stretches from the position of 
the "blob", and into the plane of the sky at a position angle of PA ^ 60°. Assuming a remnant 
age of 1000 years and the usual LMC distance, the structure has an inclination angle of about 
65° to the line of sight. The position angle is close to the symmetry axis with present and past 
activity in the visual and the X-ray wavelength ranges. We speculate that the pulsar is posi- 
tioned along this activity axis, where it has a velocity along the line of sight of a few hundred 
kms '. The "blob" is most likely a region of shock activity, as it is mainly bright in [S ii]; 
future observations of [O ii] /1/13726, 3729 would be useful to test whether the S/O abundance 
ratio is higher than average for that location in the remnant. The striking resemblance in X- 
rays between the pulsar wind nebula (PWN) of SNR 0540-69.3 and the Crab, in combination 
with our findings in this paper, suggests that the symmetry axis is part of a torus in the PWN. 
This is in agreement with the original suggestion by Gotthelf & Wang. 
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1 INTRODUCTION 

The 50.2 millisecond pulsar PSRB0540-69.3, including the sur- 
rounding supernova remnant SNR 0540-69.3, in the Large Mag- 
ellanic Cloud (LMC), is often referred to as the "Crab twin". Its 
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assumed age of about 1000 years, its pulse period, and its surround- 
ing pulsar wind nebula (PWN), are all similar to the Crab. The most 
obvious differences between the two pulsars and their surround- 
ing nebulae are that the supernova ejecta of SNR 0540-69.3 are 
oxygen-rich, and that the remnant has a shell of fast ejecta that in- 
teracts with interstellar or circumstellar g as. The latter is not yet 
confir med for the Crab (see, for example, iLundgvist & Tziamtzis 
'2012"). Whilst the Crab originates in an 8-10 M^ star (e.g.. lHestei 
2008), the prog enitor of SNR 0540-69.3 probably ha s a mass 
of about 20 Mo j William j|2olol : iLundqvist etal. II2OIIL hereafter 
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Ll l).'Serafimovich et al. '(2004') compare PSRB0540-69.3 and the 
Crab. Williams et al. (2008) review PSR B0540-69.3, its wind neb- 
ula, as well as the entire supernova remnant (SNR). 

Closer to the centre, SNR 0540-69.3 appears to be more 
asymmetric than the Crab Nebula, with much of the emission 
coming from a region a few arcseconds sou thwest of the pulsar. 
iMorse etal.1 (12006. hereafter M06; also see lKirshneret^ll989l 
and ISerafimovich et al. 1120051 herafter S05) show that this asym- 
metric appearance - as well as a red-ward asymmetry in integrated 
emission-line profiles - does not signal an overall asymmetry of 
the PWN, as the maximum velocity of emission lines towards and 
away from us are nearly equal. A detailed comparison between 
the PWNe of the two pulsars is more diffi cult due to th e large 
distance of PSRB0540-69.3 (about 5 Ikpc, 'Panagia"2005'), com- 
pared to the Crab (2±0.5kpc; see Kaplan etal. 2008, and refer- 
ences therein); the semi-major axis of the torus in the Crab only 
subtends about r.'5 at the distance of PSRB0540-69.3. Because of 
this there is still no counterpart for PSRB0540-69.3 and its PWN 
to the detailed synchronised visual and X-ray study of the Crab 
PWN by Hester et al. ( 2002). The latter authors reveal the presence 
of outward moving equatorial wisps at half the lightspeed. Many 
time-variable sub-arcsecond structures exist in the Crab PWN in 
the visual wavelen gth range i Hester et al. Il2002h . the near-infrared 
wavelength ra nge (Melatos etal. 2005), and in the X-ray regime 
dHester et al. 112002: Weisskopf et al. 2000) ; this is all discussed in 
deta il bvlHestej (|20()8|). 

Ide Luca et al. I ( 120071) are the first to report evidence for 
changes in the structure of SNR 0540-69.3 over a relatively short 
time period. They find a ten-year flux variation (1995-2005) in the 
southwest direction where the PWN emits the strongest, and sug- 
gest that they have detected a hot spot that moves at a speed of 
about four per cent of the lightspeed. They also note that the hot 
spot could be similar to a time- varying arc-like feature in the outer 
Crab Nebula, and that a pulsar jet in SNR 0540-69.3 could be di- 
rected towards the bright southwest region rathe r than perpendic- 
ular to this direction, as is earlier suggested by lOotthelf & Wand 
( l200d ). 

More recently L 11 present a study of the PWN o f PSR B0540- 
69.3; we use the same data as Ide Luca et al. I ( l2007i) . but comple- 
ment it with all available Hubble Space Telescope (HST) data. In 
particular, we include HST polarisation data from 2007, as well as 
all available X-ray data. These data confirm a remarkable activ- 
ity in the southwest region of the PWN, with the apparent "blob" 
r.'3 southwest of the pulsar, which emerges in 1999 and then 
fades until 2005-2007. We suggest that the "blob" structure may 
be a part of now fading filamentary structure that was excited in 
1999, rather than a moving entity. Particularly evident is how well 
the continuum emission traces the emission structure in [S n], and 
how low the correlation is between emission in the continuum and 
[Om]. A high intensity ratio /[Sii]//[0iii] can mean either shock ac- 
tivity ((Williams etal. 2008), or a locally high S/0-ratio. The for- 
mer scenario could, e.g., be an outflow in the PWN from the pul- 
sar, which temporally excites filaments at the "blob". Lll dis- 
card time-dependent photoionization by emission from the X-ray 
blob as the source of the enhanced /[Sii]//[0iii], and instead argue 
for slower shocks that are transmitted into PWN-embedded fila- 
ments; in a way similar to the shocks that are transmitted into the 
blobs of the circumstella r ring around SN 1987A (see, for example, 
iGroningsson et al. II2008I) . 

Unlike imaging observations, spectroscopic observations al- 
low measurements of separate emission lines and kinematics. M06, 
for example, present long-slit observations that are offset by 1" 



Table 1. The observations journal of the PSRB0540-69.3 field. 



Date 


Obs. 
Block id 


Seeing 
[arcsec] 


Ainnass 


Grism 


2010 Nov. 1 


504960 


(y.'82 


1.40 


HR-orange 


2010 Nov. 1 


504961 


(y.'84 


1.41 


HR-orange 


2010 Nov. 3 


504967 


0'.'78 


1.40 


HR-blue 


2010 Nov. 7 


504959a 


I'.'IO 


141 


HR-orange 


2010 Nov. 28 


504953 


0'.'63 


1.42 


HR-orange 


. 2010 Nov. 28 


504958 


(y.'56 


1.47 


HR-orange 


2010 Nov. 28 


504959b 


a' 60 


1.40 


HR-orange 


. 2010 Nov. 29 


504962 


0'.'62 


1.40 


HR-blue 


2010 Nov. 29 


504963 


0'.'74 


1.44 


HR-blue 


2010 Dec. 1 


504964 


(y.'75 


141 


HR-blue 


2010 Dec. 1 


504965 


(y.'75 


1.45 


HR-blue 


2010 Dec. 1 


504966a 


I'.'OO 


1.53 


HR-blue 


2010 Dec. 2 


504966b 


a'95 


1.44 


HR-blue 



west of PSRB0540-69.3, and measure various emission lines 
in the wavelength range 3335-6945A. They also present a one- 
dimensional kinematic structure along the slit (cf. their fig. 5). It 
is evident that a spatially resolved spectroscopic study in the vi- 
sual wavelength range, which completely covers the central parts of 
SNR 0540-69.3 in several emission lines, would provide new infor- 
mation. The derived velocity field could be used to create a three- 
dimensional model of the PWN, which can be used to separate and 
analyse different structures in the object. Here we present and dis- 
cuss such new observations of the central structure of SNR 0540- 
69.3. For this purpose we used the fibre-fed integral-field spectro- 
graph VIMOS at the Very Large Telescope (VLT) of the European 
Southern Observatory (ESO) in Chile. We organised the paper as 
follows. In Sections 2 and 3 we describe and discuss the observa- 
tions and the data-reduction procedure. In Section 4 we present the 
results, and follow up with a discussion in Section 5. We close the 
paper with conclusions in Section 6. 



2 OBSERVATIONS 

The object was observed in service mode during the period P86 
with the visual multi-object spectrograph ( VIMOS, iLe Fevre et al. I 
|2003) that is mounted on the Nasmyth focus of the UT3 8.2m tele- 
scope at the VLT in Chile. We used the integral field unit (IFU) with 
the smaller scale setting, 0'.'33. The HR-blue and the HR-orange 
grisms were used to observe [O in] A5007 and [S ii] AAbTlT, 673 1 . 
HR-blue (HR-orange) approximately covers the wavelength range 
4150-6200A (5250-7400A) at a resolution of i; = 2550 (i; = 2650) 
and a dispersion of 0.51 ApixeP' (0.6 A pixel"'), or 31 kms"' at 
5007A (27 km s"' at 6717A). Our object is faint, which is why we 
used the (nearly) longest exposure time possible in service-mode 
observations, 1800 s; with shorter exposures readout noise domi- 
nates faint regions with the 0'.'33 scale setting. We present our com- 
plete observations journal in Table[T] For each science exposure the 
five columns provide: the observation date, the ESO-specific ob- 
servation block id number, the seeing, the airmass, and the grism. 
The airmass and the seeing were within the ranges 1.40-1.45, and 
0'.'50-0'.'80, respectively. In this study we only used the best-seeing 
exposures that are prefixed with a bullet; we intend to combine and 
use all data in a future study. 

In addition to the science exposures (mandatory) continuum- 
lamp and arc-lamp exposures were taken immediately after each 
science exposure, as well as spectrophotometric standard-star ex- 
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Figure 1. An [O iii] /!5007-filter 40"x40" NTT/EMMI image from 1996 of 
the inner region of SNR 0540-69.3 (Lundqvist et al., in prep.). The 13"x 
13" central part illustrates region that we observed with the VIMOS IFU 
using the 0'.'33 scale setting. The pulsar (blob) position is marked by a white 
cross (circle). 



posures for each grism and night. The normalised continuum lamp 
images were used to correct for varying fibre-to-fibre transmission. 
Separate observations of the sky were considered too costly. Alto- 
gether the four detectors of the VIMOS IFU hold 40x40 = 1600 
separate fibres, where each fibre represents a spatial element, or a 
so-called spaxel, on the sky. In the 0'.'33 sampling mode each point- 
ing covers an area of I3'.'2"xI3'.'2" = 174"^ on the sky. However, the 
outermost frame of spatial elements is vignetted, which is why the 
effective area coverage is 12'.'5 x 12'.'5 = 157"^. Figure[T]shows the 
PSRB0540-69.3 region (Lundqvist etal., in preparation), as well 
as the coverage of the VIMOS IFU 0'.'33 scale setting. This smaller 
setting clearly offers the best spatial coverage and resolution. 



3 DATA REDUCTION 

Data wer e reduced using p3d (revision 15 06, which is part of ver- 
sion 2. 1. ISandin et"dni20ld |20I ll. 120121^ . VIMOS data are pro- 
vided in files of four separate detectors per field, where each de- 
tector holds 400 full spectra. The files of each detector are reduced 
separately, and are combined after flux calibration, to create an im- 
age containing the full field of all 1600 spectra. Here follows a 
description of the method. 

At first the five provided bias images of each night are com- 
bined to create a master-bias image, which is subtracted from all 
other images of the same night. In a second step, a trace mask 
is created using a continuum lamp image. Since VIMOS spectra 
overlap significantly on the CCD, cross-dispersion line profiles are 



calculated in this step, to allow a consecutive spectrum deconvolu- 
tion. Cross-dispersion line profiles are calculated every 20th pixel 
on the dispersion axis to allow an accurate extraction, they are inter- 
polated for the remaining pixels. In a third step a dispersion mask 
for the wavelength calibration is created using an arc image. There 
are few and mostly no visible arc lines blue-wards of the helium 
line at 5852 A in data of the HR-orange grism. This is no problem, 
however, since we are only interested in [Sii] /1/16717, 6731. The 
residual r between the fitted wavelength, which is calculated us- 
ing a fifth order polynomial, and the known wavelengths of the arc 
lines, is always in the range 0.02<r<0.06 A. 

Cosmic-ray hits are removed, before spectra are extracted, in 
the sc ience images using the PyCosmic algorithm JHusemann et al. I 
[2013); in the version th at is integrated into p3d, using default 
parameters. Furthermore. ISandin et al. I (120101) show that accurate 
spectrum centre positions are crucial when extracting spectra ac- 
curately, even an offset of O.I pixels affects results significantly. 
With VIMOS calculated profiles have to be ofi^set due to instru- 
ment flexure between the times of the continuum lamp exposure 
and the science exposure. New centre positions are calculated (for 
working fibres) in the science image, which is first median-filtered 
on the dispersion axis. One median value of the difference between 
all old and new centre positions is calculated and added to the pro- 
file positions. In a fourth step spectra are ext racted using the multi- 
profile deconvolution method ( Sandin et al. 12010. w hich is the p3d- 
version of the algorithm of lSharp & BirchalllCOlOl) . 

Data of all four quadrants of either grism are setup to use the 
same wavelength array. This choice allows a straightforward com- 
bination of reduced and flux-calibrated images into one spectrum 
image. The dispersion mask was then applied in a fifth step. In a 
sixth step a correction to the fibre-to-fibre sensitivity variations was 
applied to the extracted and wavelength-calibrated science spec- 
tra by dividing with the mean spectrum of an extracted and nor- 
malised continuum lamp flat-field image, excluding all spectra of 
broken and low-transmission fibres. CCD fringing effects are re- 
moved rather well in this step as the flat-field spectra are neither 
smoothed nor replaced with a polynomial, which is the usual pro- 
cedure with p3d for other instruments. During extraction data are 
changed from a CCD-based format to a row-stacked-spectra (RSS) 
format. 

In a seventh step, extracted spectra are flux calibrated using the 
reduced standard-star exposure of the respective night and grism. 
In an eighth step flux calibrated images of the four separate de- 
tectors are combined to create one RSS image with 1600 spec- 
tra. For unknown reasons, as well as insufficient information to do 
calibration even more accurately, measured flux es vary with both 
the detector and with spatial elements (also see iLagerholm et al. I 
2012). All spectra are therefore normalised to achieve the same in- 
tegrated median flux in a telluric line. We used the telluric line near 
A = 4866A {A = 6864A) with data of HR-blue (HR-orange). In 
step nine, the detector-combined images are corrected for diiferen- 
tial atmospheric refraction (DAR) as is described in Sand in etal. I 
( 2012 ). Finally, in step ten extracted RSS images are converted to 
data cubes, where the first two dimensions contain spatial informa- 
tion - East is to the left and North is up - and the third dimension 
contains the spectrum information. 



4 RESULTS 



' All papers are freely available at the p3d project web site |http ://p3d. sf .net] 
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To extract velocity structures of gases emitting in spectral lines 
that are available in our data, we focus on [S 11] AA6716, 6731 and 
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[O m] AAA959, 5007. These two multiplets produce different struc- 
tures in HST imaging observations (Lll). In particular, we find a 
spatial anti-correlation between [S ii] and [O iii] at the position of 
a "blob" about r.'3 southwest of the pulsar, with [S ii] tracing the 
optical continuum emission, as well as the surface intensity of X- 
rays. 

Potential problems with imaging observations are contamina- 
tion by point sources - including stars and the pulsar - and the LMC 
background, also the red wings of [S ii] A61i6 may blend with the 
blue wing of [S ii] /1673 1 , and only a fraction of the combined mul- 
tiplet flux enters a spectral filter. The difference in wavelength be- 
tween the two sulphur (oxygen) lines corresponds to 641 km s"' 
(2890 km s"'). In comparison, the total line width of [Om] is at 
least 3300 km s"' (M06), which means that both lines have blended 
line components. 

With our spectroscopic VIMOS data we are able to bypass 
several of these problems. At first, we note that out of 1600 spec- 
tra in one science exposure about 200 passed through broken fi- 
bres, are vignetted, or are otherwise affected by such broken and 
vignetted spectra during the DAR correction. Hence, in our anal- 
ysis we used the (approximately) 1400 remaining spectra. More- 
over, for each individual useful spectrum we analysed the full [S ii] 
and [O m] multiplets. We made linear fits to the continuum levels 
within the spectral domains of each multiplet, and then subtracted 
these continua from the respective spectra. We did not deredden the 
spectra. 

In the following we first discuss the [O in] and the [S ii] spectra 
in Sections 4. 1 and 4.2. Thereafter, we study intensity and density 
maps in Section 4.3, and velocity maps in Sections 4.4 and 4.5. 

4.1 Details of the [O ni] spectrum 

For [Om] we decomposed the two line components, us- 
ing the knowledge that the transition probability ratio 
A([Oiii]^5007)/A([Oiii]i4959) = 3. To remove the LMC 
background, we fitted the spectrally unresolved Gaussian line 
profiles in areas near the four comers of the field (avoiding spectra 
of the low-transmission fibres in the second quadrant). Thereafter, 
we created a two-dimensional fit to subtract that LMC component 
from all spectra. Simultaneously we transferred the spectrum to the 
velocity of the LMC [O iii] emission in the field, which we found 
to be vlmc = +286.5 km s"'. 

The average spectrum of all 1400 fibres in the entire field 
for [O ni] /i5007 is shown with a magenta-coloured line in the 
lower panel of Fig. |2] The full line profile fills the velocity range 
-1650 ^ I'o,,, ^ +1700 km s~', i.e., we recover the total width that 
is observed by M06; we note that our spectral resolution is 
about eight times as high as theirs (250 km s"'); their line pro- 
file could therefore be slightly broader for the same area probed. 
dLundavist & Tziamtzisll2012l discuss the importance of sufficient 
spectral resolution to constrain the faint parts of line profiles in 
the Crab SNR.) Apart from that, we cannot obtain a perfect match 
between our data and those of M06, despite their wide (2'.'0) slit 
and fairly poor (I'.' 5) seeing, which should include much of the 
PWN in the detected emission. The most obvious features are, how- 
ever, identical, i.e., the generally redshifted emission with a broad 
shoulder in the blue. What may be more pronounced in our data 
is the stronger shoulder emission in the blue at velocities vom S - 
1300 km s"' than red- wards of I'o III S + 1300 km s"'. 

Our VIMOS data allow a detailed study of spectra for smaller 
regions of the PWN. This is highlighted by the blue and the red 
spectra in the lower panel of Fig. |2] The spectrum that is drawn 
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Figure 2. Average spectra of [Sii]/ii6716, 6731 (top panel) and 
[O III] /i5007 (bottom panel) for the entire field (magenta; which includes 
about 1400 fibres out of 40x40 available fibres) and for a region that is 
centred on the "blob" that is discussed in Lll. The red (blue) line shows the 
average spectrum of 3x3 fibres (5x5 fibres). The average profiles for [S n] 
were drawn using a "master" [S n] line profile, which contains information 
of both [S n] ^6716 and [S n] i6731. 



in red colour is the average spectrum of 3 x 3 fibres (correspond- 
ing to a l"x 1" area), which encapsulates the "blob" in Lll. The 
spectrum that is drawn in blue colour is the average spectrum of 
5x5 fibres, which are also centred on the "blob". In general, the 
average spectra of the "blob" are much stronger than the average 
spectrum of the entire field, and there is nearly no emission near 
v'oiii = +100kms"'. This low emission is flanked by peaks at about 
vo,„=-250kms"' and I'om = +450kms"', as well as a second peak 
at vo III = +900 km s"'. Compared to the 1400-spectrum average line 
profile, the "blob" spectrum is centred to the red, with a possibly 
higher red shoulder at Vq,,, > + 1300 km s"'. The blue shoulder at 
velocities Vo,„ < -1200 ±300 km s"' is very similar to that of the 
spectrum of the entire field. This suggests that the blue shoulder is 
a general feature that exists across the entire field. We discuss this 
further in Section 4.4. 



4.2 Details of the combined [S n] spectrum 

The LMC background is lower at the emission lines of [S n] than 
at [Om], but a complication with recovering a clean [S n] line pro- 
file is the line ratio C[s,,j =/([Sii] ^6731)//([Sii] .16716), which de- 
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pends on the electron density. 805 identify a region wiiere the two 
line components do not overlap, and use this region to derive a line 
ratio and then the electron density ^e = (1 - 5) x 10^ cm"'. Here 
we decomposed the two line components for all spectra where the 
signal-to-noise is high enough. Next we describe the procedure. 

After we removed the continuum emission in the same way as 
for [O III], we fitted the line spectrum (before decomposition) using 
a Savitzky-Golay filter to mimic noise influence. We then swept 
from the blue side using [S n] /16716 as a template, and from the 
red side using [S ii] A6731 as a template. From cisuj we constructed 
a "master" [S ii] line profile f^: 

/a = w/ii +{l -w)fA2, W = , 

1 + C[S „] 

where /ji is the forward sweep for [S ii] /i6716 and /o the back- 
ward sweep for [S ii] /i673 1 . For low-density gas dominating the 
spectrum, [Sii]/i6716 dominates the "master" [Sii] line profile, 
and for high-density gas [S ii] /16731 dominates. We calculated C[s„] 
to give the best correlation between the blue-ward sweep using 
[S ii] /16716 as template, and the red- ward sweep using [S ii] /16731 
as template. The average correlation coefficient between those 
functions for the central 34x34 fibres is 0.90. For the 5x5 fibres for 
the "blob" in Fig.[2l the average correlation coefficient is 0.95. 

To remove the LMC background, we used a similar procedure 
as with [Om]. The velocity of the LMC [Sn] emission, was also 
in this case found to be Vlmc = +286.5 km s"'. The average of all 
1400 spectra is shown in the upper panel of Fig. |2](magenta). The 
full [S ii] line profile is not as wide as the corresponding line profile 
of [Oiii]/I5007, -1200^vis„]s;i200kms-'. Compared to [Om], 
the [S n] line profile is more heavily biased towards redder wave- 
lengths, with a pronounced peak at V[s„] ~ +300 km s"'. As with 
[Om], there could be a shoulder in the blue, but this only extends 
out to V[s II] =- 1 200 km s" ' . 

The blue and the red spectra in the upper panel of Fig. |2]high- 
light the area around the "blob". The difference in emission per 
spectrum between the "blob" region and the entire field is even 
more pronounced for [S n] than for [O m]. There is also here nearly 
no emission for ^ Vfs,,] ^ 100 km s"', and there are peaks at 
V[sii] ~-300kms"' and V[Sii] ~ +400 kms"'; although, unlike for 
[O III] there is no peak at V[s „] ~ +900 km s" ' . 



4.3 Intensity and density maps of both [O m] and [S n] 

We integrated individual spectra of [Om] and [Sii] in velocity 
to generate projected intensity maps, see the upper two panels of 
Fig-El To highlight differences in the surface-intensity distribution 
of [Om] and [Sii], we show the intensity ratio /[Sii]//[Oiiij in the 
lower right-hand side panel; again, note that the spectra were not 
extinction corrected. We also show an electron-density (A'e) map in 
the lower left-hand side panel, which is derived from the decompo- 
sition of the two line components of [S n]. To calculate the electron 
density, we used o ur multilevel model atom that is described in 
lMaranetal.1 ( l2000h . 

A few regions stand out in a comparison of the maps of 
[Om] and [Sii]. To the northeast, there is a region where both 
/[s II] //[o III] and N^ are low. A similar region exists immediately west 
of the "blob", where A'e is higher. The intensity ratio /[Sii]//[Oiii] 
is high for the "blob" itself, and in the northwest region of the 
PWN. Note the low density of the "blob". Moreover, high densities 
are seen in the southeast, but neither [O iii] nor [S ii] are particu- 
larly bright there. In general, we do not recover the high value of 
A'e = (1 - 5) X 10' cm"' of S05. Instead, the electron density in the 



[S ii] filaments is < 10' cm"', the highest value is about 750cm"' 
< I" west of the blob and 1-3" east of the blob. We caution, how- 
ever, that the maps in Fig.|3]were smoothed over two pixels, which 
results in smeared out density peaks. In unsmoothed data there are 
single spatial elements where A'e - 2x 10' cm"', and in this sense 
the high value of S05 is recovered; their value is measured at one 
position along their slit where the two components of [S ii] do not 
blend. 

In a comparison with the wavelet-filtered HST images in 
fig. 13 of LI 1, we recover the general structures of the HST [S ii] 
and [O m] images well. The difference in epoch between the HST 
and the VLT data sets is about eleven years. A peak in /[Sii]//[Oiii] - 
0.4 occurs at the "blob" position in both data sets. Another peak 
with the same value on the ratio is seen in the northwest region, 
where the remnant bulges out in [S n] (cf. Fig.[3l(; an enhanced ra- 
tio is also present at that position in the HST intensity-ratio plot. 



4.4 Intensity maps at separate velocities 

Our spectroscopic data contain more information than the projected 
images in Fig.[3]reveal. We show intensity maps of separate veloc- 
ity intervals for [O m] in Fig. [4] and for [S n] in Fig. [5] These fig- 
ures are smoothed in a convolution using a 3 x 3-pixel kernel. We 
only show the central 32x32 spatial elements (10'.'56x ^.'56 or 
2.61 x2. 61 pc^, assuming LMC is 51 kpc away, cf. Section lSTt . 

We begin with gas that moves the fastest towards us - the 
three panels in the uppermost row of Fig.|4]for [O m]. There is an 
extended area of emission, with some enhancement, in the north- 
east. This emission strengthens our hypothesis in Section 4. 1 that 
the blue shoulder of the 1400-spectrum and the "blob" spectrum 
in Fig. [2] possibly is an extended fast shell-like structure in [Om] 
on the near side of the remnant, as seen from the Earth. A corre- 
sponding emission, or "haze", is visible at high positive velocities 
(Fig. [8}, whilst there is no clear red shoulder in Fig. |2] The "haze" 
is thus weaker on the red side than on the blue side. 

The "blob" emerges in the velocity interval -800 ^ 
V ^-300kms"', in both [Om] and [Sii]. A pronounced region 
southeast of the "blob" is also visible. For -300 ^ v ^-100 kms"', 
the "blob" becomes encircled by a shell that breaks up at higher 
(redder) velocities; the "blob" appears to move at a lower veloc- 
ity. A point-like structure re-emerges at the "blob" position for 
400 ^ V ^ 550 km s"'; this explains the line profiles in Fig.|2] The 
main difference between [O m] and [S n] is, furthermore, that [S ii] 
is also very pronounced for 250 ^ i' ^ 400km s"', whereas [Om] 
peaks immediately west of the "blob". Figure|4]shows [Om] emis- 
sion from the "blob" for 850 ^ v ^ 1400 km s"' , which confirms the 
high- velocity peak in [Om] in Fig. [2] This structure also appears to 
be connected to a region northeast of the "blob". For the "blob" it 
seems as if it could be a separate entity for -800 ^ i' ^ 700 km s"' . 
The high value of /[Sii]//[Oiii] for the "blob" in Fig. [3]is likely due 
to strong [S ii] emission for 250 ^ v ^ 400 km s" ' , where [O m] is 
fainter. 

Earlier in this section we highlighted additional regions. The 
region to the northeast, with a low value on both /[Sii]//[0iii] and A'e, 
appears to be dominated by [Om] for 100 ^ v ^ 250 km s"'. The 
(high-density) region immediately west of the blob, which shows a 
low value on /[Sii]//[Oiii), emits strongly in [Om] at many velocities 
in the range -300 ^ v ^ + 1 100 km s"' . In the same region, the west 
streak of high electron density and high /isii]//[Oiii] in Fig. |3]- that 
extends from north to south between the protrusions in [S ii] to the 
northeast and the southeast - is marked mainly for -100 ^ v ^ 
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Figure 3. For the innermost 1 1'.'22 X 1 1'.'22 (2.77x2.77 pc^) area of SNR 0540-69.3, tlie flux of [S n] ([O iii]) is shown in the upper left-hand side (right-hand 
side) panel. The intensity unit is intensity per (X'33 X 0'.'33 pixel. The images show no continuum sources, since we subtracted the continuum for each spatial 
element individually. Data were smoothed in a convolution with a 2x2-pixel kernel. The lower left-hand side (light-hand side) panel shows the electron density 
^e (^[Sii]//[Oiii])- Noith is up and east to the left. Each spatial element is 0'.'33xO'.'33. Spatial elements around the border are differently aifected by effects 
of the DAR correction. The spatial elements in the northeast region show emission that passed through low-transmission fibres in the second quadrant. The 
contours of the [S 11] flux are, moreover, reproduced in all panels. The position and size of the "blob" (LI 1) is indicated by a magenta-coloured circle and the 
pulsar position is marked by a white X. 



250 km s ' . Despite a high density, this feature could partly be due 
to poorly subtracted LMC emission. 



4.5 A closer study of the velocity structures 

In Fig. [6] we present projected intensities of the three-dimensional 
data cube of [Out] and [Sn]. The intensities were summed along 
the line of sight (the .i;-axis; the y-axis) in the upper (middle; lower) 
panel of the figure, where the x-axis (velocity-axis; velocity-axis) 
extends rightwards from East to West (approaching-to-receding; 
approaching-to-receding) and the v-axis (v-axis; j-axis) upwards 
from South to North (South to North; East to West). The two up- 
per panels show the lines as they appear on the sky. To guide the 



eye we put a raster on top of each image. The "blob" is centred 
at {x,y) X (22.1, 13.7). The images were not smoothed, unlike the 
data in Fig. [3] To see details in the central part of the remnant, we 
focused on the velocity interval -1200 ^v^ 1200 km s"'. The ab- 
solute intensity-scale panels in Fig.|6]complement those in Figs.|4] 
and [5] since the panels in the latter two figures use different inten- 
sity cuts to highlight various features in the individual panels. 

In the second row of Fig. |6l we clearly see three distinct ve- 
locity features: one near the LMC rest velocity, another broader for 
200^ v^ 500km s"', and a third one near v = 800kms '. The first 
feature is located to the northern part of the data cube and is most 
likely influenced by the LMC background. Focusing on the struc- 
tures between (5,5)^ (x, y) ^ (30, 30) (the region is indicated with 
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Figure 4. Intensity maps of [O ni] at separate and specified velocity intervals for the innermost 10'.'56x 10'.'56 (2.61 x2.61 pc^) area of SNR 0540-69.3, each 
spatial element is 0'.'33xO'.'33. The position and size of the "blob" (Lll) is indicated by a magenta-coloured circle and the pulsar position is marked by an X. 
Colours are individual to each panel to emphasize structures in the selected velocity range. 
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Figure 5. Intensity maps of [S ii] at separate and specified velocity intervals for the innermost 10'.'56x 10'.'56 (2.61 x2. 61 pc^) area of SNR 0540-69.3, each 
spatial element is 0'.'33xO'.'33. For more details see the caption of Fig.|4] 



a large white box in both top-row panels), it is clear that some of 
the emission about zero velocity at higher x values may not belong 
to the remnant. The situation is different for 5 ^ y ^ 15, where 
there is a blueshifted feature in [Out] at about v =-200 km s"', 
which also continues out to v»-500kms"' (indicated with a white 
box in the middle left-hand side panel). For both [Oiii] and [Sii], 
there is enhanced emission about the y-value for the "blob", i.e., 
(18, 1 1) ^ (x,y) < (25, 17) (the "blob" is indicated with a white box 
with rounded comers in both top-row panels; this is more easily 
seen in the on-line version). Furthermore, in the third row of Fig.|6] 
the "blob" region clearly sticks out; [S n] emission from the "blob" 
region begins at about v =-500 km s' and ends with bright emis- 
sion at about v = 400 km s"' (indicated with a white box). This is 
also seen in the velocity plots Figs. |4] and |5] where the "blob" is 
clearly seen for -800^v^-300kms"'. 

The second velocity component in the second row of panels, 
200 < V s; 500 km S-' , subtends the interval (7, 8) < (x, y) ^ (27, 27) 



(using the third-row panels; the region falls just inside the large 
white box in the top panel). In a comparison with Figs. |4] and |5] we 
clearly see that there are mainly two emitting regions for 200 ^ v ^ 
500 km s" ' . In [O m] the two regions are (20, 1 2) s^ (x, y) ^ (25, 1 8) 
and (7, 15) ^ {x,y) ^ (12, 25) (these two regions are indicated with 
dash-dotted boxes). In [Sii], the former region dominates; this is 
clearly seen in the lower right-hand side panel of Fig. |6] Finally, the 
third component, at about 800 km s"', is fainter and mainly emits 
in [On:]; it occupies the interval (7, 12) $; {x,y) ^ (19,23) (indi- 
cated with a dotted box). This component is not the strongest at the 
"blob" position, but in the eastern region of the PWN. 

An interesting feature of Fig. |6] is the diagonal ubiquitous 
structure that extends in [Oiii] from x = 20 to x =: 12, between 
V X Okms"' and v » 800 kms"', respectively, for the three major 
velocity features in the lower left-hand side panel (indicated with 
a white solid line). In the middle left-hand side panel, the diagonal 
is less obvious, but if we neglect the northern part of the velocity 



© 2013 RAS, MNRAS OOQ.fTHTs] 



The three-dimensional structure ofSNR 0540-69.3 9 



[om] 



4t-l4 



J,^U 




4<!-L4 
3. Jr. 14 
k N 

2u-M 
SJf-t4 

a 



ic I.-: SI :.■ 



[om] 






If- I*- 




aild He- lb 

3Ja-l« 

3e.-t6 



lODft .SO!) n 5nti 

Velodty (km ^ £.1 



IflOO 



[o mj 



■ 


^^^^ 






40 


^^^^^H| 






■ 




^^^H 




.^! 


■ 


^H 






■ 


;Ki 




H 






■ 








It 


11 






H 








3(1 






I 


..S 


^^■1 






10 


■ 


^D 


B^W ._ 




■ 


^1 


^ y 





4b^1« — 













3e-J& 


— 




I<;-1* 


r.<^!«i 




1 


k If. 
J:-lli 




1 


|(3c-lt. 

Ml? 
I> 


f 



Veli3Qbr(km,'£l 



[sni 



st-i; 
VIS 




Ij Z\'t 1^ 3i:- 5.'- *J 



[Sill 



lu-.t* ■ 
Scl-J- . 



i ^^^^^H 



■4 



l(KO .'iDO n SCO 

V«k);l!i"'yjn/s) 



1000 



[SII] 

(in1rTis.it>' in crip /sytmJ / (km .' ^1) 



IJirfS- 



li-(S 

, 4H!-I'f 

it IT 

■ 




LSr-lft 

>:^ I.4E.IS. 
l.-lr-JB- 

ilT 
' ii 



jm SM 



Figure 6. Integrated data of three-dimensional cubes in [O m] (left-liand side column) and [S n] (right-hand side column) as they appear on the sky (upper row). 
Increasing values on the x-axis (y-axis) go from East to West (South to North). The second (third) row is obtained by summing values along the East-West 
(South-North) direction. White lines and boxes indicate regions that are discussed in the text. A raster has been put on the panels for oiientation; one cell is 
62.5kms"' wide on the velocity axes. The "blob" is centred on (x,y)K:(22.l, 13.7). 
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structure at the LMC rest velocity, we can draw a diagonal from 
y =^ 14 to y =: 18 between v--200kms"' and v == 800 km s"' (in- 
dicated with a white solid line). The same diagonal structures also 
appear for [S ii]. The "blob" is projected at the low-velocity ends of 
those two diagonals. Also note the symmetric features about those 
endpoints in velocity, for v ^ ±400 km s"' . 



5 DISCUSSION 

In the following we model three-dimensional structures of 
SNR 0540-69.3 in Section 5.1, to interpret the new data. There- 
after, in Section 5.2 we compare our observations with previous 
studies. 



5.1 A closer study of the three-dimensional structures 

In the previous section we found hints for diagonal structures, 
which originate in major features in the remnant. We speculate that 
such diagonal structures are part of bulk motions of the expand- 
ing gas. To explore this possibility further we used information of 
these ejecta to build a three-dimensional model of the remnant, as- 
suming a linear expansion. Additional assumpti ons w ere an age 
of 1000 years and a distance of 51 kpc ( Panagiall2005h . This im- 
plies that each 0'.'33xO'.'33 spatial element on the sky corresponds 
to 0.082 X 0.082 pc^ and that 1000 km s"' along the line of sight is 

1.02 pc. Because the seeing was about 0'.'6, the spatial resolution is 
approximately 0.15 pc. With a velocity resolution of 31 kms"'(cf. 
Section 2), the resolution along the line of sight is 0.032 pc, i.e. 
about 4.6 times as good as the spatial resolution, and only a factor 
2.8 worse than the best spatial resolution of /f5'r/WFPC2, 0.01 1 pc, 
obtained with the Planetary Camera chip used in Lll. We defined 
the X and the y axes as in Fig. [6] and added a z axis that is oriented 
towards the observer, along the line of sight. We emphasise that 
the real structure along the z-axis - as derived from motions along 
the line of sight - may not be the same as the one we derive. For 
example, if gas emitting [O in] or [S ii] was accelerated throughout 
the evolution or suddenly, as discussed in Lll, structures may be 
exaggerated. However, if we focus on the larger features of the rem- 
nant, we could be close to the real structure. As another example, 
filaments in the Crab expand nearly linearly, and a linear expan - 
sion assumption is correct to within about 10% ( lTrimblelll968l) . 
The main uncertainty in a transformation from radial velocity to a 
physical scale for SNR 0540-69.3 is probably the uncertainty of 
the remnant age. 

With our derived properties of the diagonal structures in Sec- 
tion 4.5 and the above assumptions about the transformation to 
physical scales, the diagonals correspond to an angle of about 65° 
relative to the line of sight. We therefore expect an overall structure 
to extend northeast from about zero velocity, and into the plane of 
the sky at an angle of 65° relative to the line of sight. The structure 
crosses the xy plane roughly at the position of the "blob", and the 
position angle relative to the sky is PA^ 60° . These angles should 
only be considered as guidelines to the orientation of the overall 
remnant structure, as the three main velocity structures in Fig. |6] 
(middle and bottom panels) are extended in both the x and the y 
directions; this leads to some uncertainty in the exact orientation of 
the diagonal structure. 

To better interpret the remnant structure, we used the visu- 
alisation tool MayaVi. We created interactive three-dimensional 
maps with this tool, using Python in an approach that follows 



IVogt & Dopitl ( I2OIOI) for S NR IE 0102.2-7219 in the Small Mag- 
ellanic Cloud and IVogt & D opita (200) for SNR N132D in the 
LMC. The data cube was read in a four-dimensional array and 
was visualised point-by-point as partly transparent spheres; data 
of columns that contain spectra of broken fibres were removed in 
advance. The flux of each point is represented by both size and 
colour. A brighter colour and a larger size of a sphere correspond to 
a higher flux at a specific point. Foreground spheres are somewhat 
transparent to allow background spheres to shine through. We show 
[O III] and [S 11] structures for three projections in Fig. |7] To dis- 
criminate between approaching and receding gas, we colour-coded 
the gas differently, with a blue colour bar for approaching gas - i.e. 
positive z - and a red colour bar for receding gas. Zero velocity 
corresponds to the local LMC rest velocity. We emphasise that the 
visualisation is mainly meant for highlighting structures, it must 
not be used in quantitative measurements. 

We used MayaVi in three steps to produce the panels in Fig. [7] 
In the upper panels of Fig. |7] we show scenes that are similar to 
those in the upper panels of Figs. |3] and |6] Intensities for negative 
and positive z-values were produced separately. We also visualise 
the fainter flux around the central part separately, see Fig. [8j this 
component has similar colour coding. The scenes in the lower two 
panels of Fig.|6]appear smaller in [S n] than in [On:]. This is merely 
an effect of the shorter velocity interval where we could measure 
line emission in [S 11]. 

The right-hand side middle panel in Fig. |7] shows an overall 
diagonal structure that extends in the plane of the sky from the 
"blob" position at positive x-values at large inclination angle with 
the z-axis. This recovers the diagonal with 65° relative to the line 
of sight that we inferred from Fig. |6](cf. above). The structure ap- 
pears to form two ellipse-like structures in [O iii], one at about zero 
velocity, and the other towards the lower left comer in the middle 
left-hand side panel. The two ellipses are plausibly ring-like fea- 
tures at the ends of a bipolar structure. According to the velocities 
in Fig.|6l the upper right-hand side structure lies further in along the 
j-axis than the lower left-hand side structure. The [S n] structure is 
more complex; the v » ±400 km s"' structures, which we mention at 
the end of in Section 4.5, are clearly visible as the global diagonal 
structure and differ from the [O m] structures. In the bottom panels 
of Fig. [7] we view the possibly bipolar structure closer to its axis as 
viewed in the plane of the sky with the x-axis oriented inwards from 
East to West; different parts overlap more, and it is more difficult to 
disentangle individual features. 

In Fig.[8]we note abundant high- velocity gas in [Om] that fills 
the boxes in all panels; in particular on the near side of the box, as 
seen from the Earth (blue colour coding). This corresponds to the 
high- velocity shoulder that is seen in the line profiles in Fig.|2] 

5.2 The structure of the pulsar wind nebula 

The large structure that we identified in the previous subsection and 
in Figs.|6]and|7]has a position angle relative to the sky of PA=:60°. 
There is certainly some uncertainty to this position angle that might 
be high enough to accommodate the symmetry axis that is dis- 
cussed in LI 1. In this previous study we find that many parameters 
change along this symmetry axis, and we argue for present and pos- 
sible past shock activity along it; the symmetry-axis position angle 
isPA=45°-50°. 

The "blob" region stands out in Figs. |6] and [7] An explanation 
to its existence could be a shock that propagates at v » ±400 km s" ' 
along the line of sight. For normal cosmic abundances, the main 
problem with this inteipretation is that cooling and recombination 
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Figure 7. Three-dimensional scenes of the central region of SNR 0540-69.3 that were created using the visualisation program MayaVi. The flux at each point 
is represented by both size and colour of partly transparent spheres. Zero velocity is at the local LMC rest velocity (i'lmc = +286.5 km s"'). Blue colours 
indicate approaching gas at negative velocities - i.e. positive z - and red colours indicate receding gas at positive velocities. Both colour bars are shown on 
the left-hand side of respective top panel. The three panels in the left-hand side (light-hand side) column show [O in] ([S n]) from difl'erent angles. The upper 
panels show the data cube as it appears on the sky; the jc-axis (y-axis) extends rightwards from East to West (upwards from South to North), and the z-axis 
towards the Earth for decreasing velocities. The middle panels show the data cube as viewed in the plane of the sky from South to North inwards on the 
y-axis; the .v-axis (z-axis) extends lightwards from East to West (upwards from longer wavelengths to shorter wavelengths). The lower panels show the data 
cube as viewed in the plane of the sky from East to West inwards on the jc-axis; the y-axis (z-axis) extends rightwards from North to South (upwards from 
longer wavelengths to shorter wavelengths). The pulsar is identified with a black sphere in all panels; it was placed at the LMC rest velocity, using the exact 
optical coordinates of PSRB0540-69.3. The vertical black line strikes through the pulsar and is parallel to the line of sight. The true position of the pulsar is 
somewhere along this line. The "blob" is located to the region that is void of "blue" emission in [O in], just southwest of the pulsar (top left-hand side panel), 
and overlaps with the region of "blue" emission in [S n] (top right-hand side panel); also compare with the projected intensity plots in Fig. [3] The axes of the 
two upper panels are contracted to leave space for the intensity bars. 
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Figure 8. Three-dimensional scenes of the [O in] faint emission in the central region of SNR 0540-69.3 that was visualised and presented in the same way as 
in Fig.|7] The panels represent different orientations of the data cube. The upper left-hand side panel shows the data cube as it appears on the sky; the ,t-axis 
(y-axis) extends rightwards from East to West (upwards from South to North), and the z-axis towards the Earth for decreasing velocities. The upper right-hand 
side panel highlights the structure of faint emission. The lower left-hand side panel shows the data cube as viewed in the plane of the sky from South to 
North inwards on the y-axis; the A-axis (z-axis) extends rightwards from East to West (upwards from longer wavelengths to shorter wavelengths). The lower 
right-hand side panel shows the data cube as viewed in the plane of the sky from East to West inwards on the x-axis; the y-axis (z-axis) extends rightwards 
from North to South (upwards from longer wavelengths to shorter wavelengths). The pulsar is identified with a black sphere in all panels. The vertical black 
line strikes through the pulsar and is parallel to the line of sight. 



times are long for densities that are typical of this object. The post- 
shock temperature would be about 2 x lO^K, and coohng times 
would be comparable to the age of the remnant for a post-shock 
electron density of lO^cm"'. The situation is, however, different 
for SNR 0540-69.3 that has a high meta l abundance. We have 
used our plasma code that is discussed in ISorokina et al. I ( 120041) 
and Mattila etal. (2008) to estimate cooling times for a range of 
ab undances. In particular , with the relative abundances suggested 
by Iwilliams et al. I ( l2008h for the oxygen-rich filaments - i.e., O: 
Ne: Mg: Si: S: Ar: Ca: Fe = 1: 0.2: 0.1: 0.1: 0.1: 0.1: 0.1: 0.1 by 
number - the cooling time scale is less than a few years for a post- 
shock electron density of 10 cm"^^ and temperatures corresponding 
to the shock velocities we tested, i.e., between 50-500 km s"'. We 
can even allow H and He to dominate, with H: He: O = 10: 2: 
1 - with all other metal abundances as in the oxygen-rich case - 
and still have a cooling time scale of less than 10 years for a post- 
shock electron density of 10 cm"'. It therefore seems plausible that 
a radiative shock that propagates at va ±400 km s"' along the line 
of sight could exist at the "blob" position, but we cannot rule out 



slower shocks than this; if they are driven into already expanding 
structures, such as the ring-like features that are mainly seen in 
[Oiii]. Observations that sample [On] emission at high spatial res- 
olution would provide useful data to constrain the shock velocity. 
Unfortunately, VIMOS/IFU does not contain [O ii] /1/13726, 3729; 
previous observations show that these lines are strong in the central 
region of SNR 0540-69.3 ( Kirshner et al. 1989; S05). 

To explain shocks that propagate at velocities of up to 
400 km s"' along the line of sight, we need to use similar argu- 
ments as Lll; i.e., some event occuiTed at the "blob" position, 
which caused shocks to be driven into filaments there. The high 
/is,i]//[0iii] ratio is not due to high local sulphur overabundance, if 
[S ii] is concentrated to this location. The reason to the local re- 
lease of energy at this location still needs to be explained, although 
energy input from pulsar activity seems plausible. 

The faint [Oiii] emission in Fig. [8] seems stronger on the near 
side of remnant than on the far side. Comparing the lower right- 
hand side panel of Fig. [8] with the lower left-hand side panel of 
Fig.|7]it appears as if the faint [Oiii] emission is particularly strong 
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at large v- values on the near side and small v- values on the far side. 
This could signal some asymmetry in the presumably photoionised 
ejecta in elongated structures along the jet axis, as in the Crab, or 
in ejecta beyond the PWN/ejecta interface. 



5.2.1 A discussion of the pulsar-activity axis 

Here we discuss a possible common origin to the structure we 
find here and the structures of LU. As we mentioned in the pre- 
vious section, we assumed that the ring-like structures - that are 
seen in Fig. |7] - expand linearly. The position of the pulsar is 
(x, y ) - (19, 17), which co rresponds to the position projected on the 
sky JSerafimovich etal. 112004) . We are tempted to suggest that the 
pulsar is positioned close to an axis that is put between the centres 
of the two ring-like features. This places the pulsar at a velocity 
along the line of sight of v == -1-400 km s"'. A 3cr upper limit o n 
the proper motion on the sky is 290 km s"' jde Luca et al. 1120070 . 
which implies that the true pulsar motion is 400-500 km s"', di- 
rected away from us at an angle that is within =: 36° of the line of 
sight. An argument against the pulsar sitting on the far (red) side of 
the nebula is the structure of faint [O in] emission (Fig. [8](, which 
appears fairly centred on the pulsar at the LMC rest velocity. 

In Fig. [7] [S n] appears to fill the structure more along the 
axis between the ring-like features than [O m]. If this axis is close 
to a pulsar-activity axis, such as a pulsar jet, or a part of a torus 
(see below), there is the possibility that some of the [S ii]-emitting 
structures could be more affected by the pulsar than structures that 
primarily emit [Om]. The "blob" could be such a [Sii] -emitting 
structure that was shocked by the pulsar activity. In such a sce- 
nario - where past and present shocks are allowed along the pulsar- 
activity axis (LI 1) - it is occasionally necessary to bring in gas to- 
wards this axis, or the pulsar in put may vary along the activity axis. 
iBucciantini et al. I ( 12004) and iStone & Gardinen ( 120070 discuss pro- 
cesses that create and pull filaments into the PWN, but it is unclear 
how efficient such a filamentation is, and how late in the PWN evo- 
lution filaments could develop to remain unsuppressed. If it is the 
pulsar activity that varies, filaments may not have to move towards 
the axis. In this case activities in the PWN along the pulsar-activity 
axis vary instead. 



5.2.2 A comparison with the Crab using a torus model 

In a comparison with the Crab PWN JMori etal. I ( |2004|) show that 
the outermost parts of the Crab PWN X-ray torus reach beyond 
50" away from the pulsar. The outermost regions of the PWN torus 
of PSRB0540-69.3 (if similar to that of the Crab) would reside at 

> 2" from PSRB0540-69.3, if we assu me a distance of 2 kpc to 
the Crab and 50 kpc to SNR 0540-69.3. iMorietal.l ( |2004|) iden- 
tify the entire torus in the Crab PWN to be external to the PWN 
termination shock, so the brightest parts of the X-ray torus should 
be in the post-shock flow. In particular, one can identify a bright 
feature at about 30" - which corresponds to about V!2 at the dis- 
tance of SNR 0540-69.3- southwest of the Crab pulsar, which also 
shows the hardest X-ray spectmm in the PWN (the photon index 

< 1.8, cf. fig 3 in Mori etal. 2004). This is reminiscent of the 
situation in the PWN of PSRB0540~69.3. In LI 1 we find that the 
region close to the "blob", r.'3-r.'4 southwest of the pulsar, shows 
the hardest X-ray spectrum with a photon index of about 1.65. In 
the Crab, this region is on the far side of the torus, and is there- 
fore not Doppler boosted. If we identify the "blob" region in the 
PWN of PSRB0540-69.3 in the same way - that is, as an active 



region in the PWN torus - we can speculate that the PWN torus 
in SNR 0540-69.3 could be Doppler boosted on the south-western 
side. The X-ray emission is namely brighter on that side compared 
to the north-eastern side. X-ray emission could also be enhanced 
due to interaction with supernova ejecta. 

If we identify the symmetry axis in the PWN of SNR 0540- 
69.3 as a torus, we note that the difference between it and the Crab 
is that the optical line emission follows the X-ray emission closely 
in SNR 0540-69.3. There is no such strong cor relation in the Crab 
(see, for example, fig. 11 in lCadez etal. 1120040 . This observation, 
and the in general higher temperatures that are derived from [O in] 
in the PWN of SNR 0540-69.3 compared to the Crab, favour a sce- 
nario where optical emission to a larger extent is driven by shocks 
rather than by photoionization in the PWN of SNR 0540-69.3. The 
true ejecta distribution in the latter object may therefore not be cor- 
rectly described by the structure as visualised in Figs.|7]and[8l since 
regions that are not excited by shocks will be fainter. The wide 
wavelength-range observational study of ICharlebois etal. I ( 120101) 
of the entire region of the Crab can provide a clue to the over- 
all structure also for the PWN of SNR 0540-69.3. Their three- 
dimensional movie shows bright loops that are visible preferen- 
tially in the torus plane, whereas structures in the jet directions are 
less structured and fainter If the distance between the torus and 
the loops were shorter, a situation may arise s uch as that in the 
PWN of SNR 0540-69.3. The discussion in Wilhams et al. I J2008h 
that oxygen-rich ejecta in SNR 0540-69.3 may have high-density 
structures that lag behind the general PWN/ejecta interface fits such 
a scenario. 



5.2.3 A comparison with Cas A, SN 1987 A, and 3D models 

If the ejecta structures we observe in [S ii] and [O m] correspond 
to the real matter distribution - and not just an effect of excita- 
tion by, for example, a PWN torus - this could lead to clues about 
the supernova explosion. A parallel case with apparent asymmetric 
ejecta structures is SN 1987A, where observed inner ejecta lie in a 
plane that is close to the equatorial plane of the circumstellar ring 
(Lar sson etal. 1120130 . Here structures of the metal-rich ejecta can 
be probed out to about 3500 km s"' , as a consequence of still active 
radioactive decay and ongoing interaction with the circumstellar 
ring. This higher velocity limit is twice as high as what we achieve 
with our VIMOS data for SNR 0540-69.3. A direct comparison 
with SN 1987A is therefore not straightforward. 

Another young well-observed c o re-collapse supernova 
remnant is Cas A dFesenetalTI 120061; iDeLanev etal.l |2010|; 
iHwang & Laminall2012l) . A comparison with SNR 0540-69.3 is 
also in t his case complicate d, as the progenitor of Cas A was a type 
lib SN dKrause et al. II 198Sl) with a n ejecta mass of perhaps as low 
as 3 Mo dHwang & Laminall2012l) . The ejecta of SNR 0540-69.3 
are more massive, probably more similar to the ejecta mass of 
SN 1987A (Lll). Due to the low ejecta mass, in combination 
with pre-supemova mass loss, the observed ejecta structures 
in Cas A are affected by the reverse shock; the reverse shock 
has now reached freely expanding ejecta at about 5000 km s"' 
dHwang & Laming 120120 . A massive asymmetric circumstellar 
shell is needed if a reverse shock were to be observationally impor- 
tant also in SNR 0540-69.3. If we envision a circumstellar disk, 
as in SN 1987A, this may still leave an imprint if it was massive 
enough. A s cenario with a circum stellar ring in SNR 0540-69.3 is 
discussed bv lCaraveo etal. l ( 1998), but it is not clear why it would 
be oxygen-rich, which is needed to explain the line emission in 
SNR 0540-69.3 (see, for example, fig. 6 in M06). 
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The most advanced three-dimensional supernova simulation 
that is available fo r the full evolu t ion fro m explosion until shock 
breakout is that of lHammer etal. I ( 120101). It is based on a 15 Mq 
blue supergiant progenitor model of IWooslev et al. I (119881 see also 
iKifonidis et al. 2003). Such a low-mass compact progenitor is most 
likely not a perfect model for the progenitor of SNR 0540-69.3. In 
addition, the three-dimensional simulations do not include radioac- 
tive heating, which underestimates the volume of the radioactive 
elements. For SNR 0540-69.3 there is also the interaction between 
ejecta and the PWN, which alters the original structure. Neverthe- 
less, we can probably draw some general conclusions on the three- 
dimensional structure of elemental abundances. The simulations 
clearly show how high-Z elements such as nickel (^*Ni) - as well 
as other iron-group elements and silicon - are mixed far out into 
the envelope, and are actually more mixed than oxygen, neon, and 
magnesium, on average. Looking in more detail, there are, how- 
ever, large diiferences in the three-dimensional distribution of the 
elements. Most spectacular is the high-Z element distribution, with 
rather few long fingers protruding into the H/He envelope. This nat- 
urally explains the necessity of artificial mixing to large velocities 
of ^^Ni in models of SN 1987A to explain its light curve (e.g., 
iBlinnikov et al. II2000I) . Also, hydrogen is mixed in to low veloc- 
ities (^ 3000 km s"'), which can explain why we see hydrogen in 
SNR 0540-69.3 (S05). 

We expect to take the study presented here further in a future 
paper. With a full analysis of all spectra in Table 1 , we expect to 
disentangle the three-dimensional structure of SNR 0540-69.3 at 
a higher signal-to-noise ratio, which will enable us to obtain a full 
three-dimensional map of the [S ii] density, and not just projected 
densities. We will also include more elements, such as argon, to 
constrain the elemental distribution. This will allow us to make a 
more detailed comparison between the derived three-dimensional 
structure of the inner eje c ta of SNR 0540-69.3 with models such 
as that of lHammer et al. I ( 1201(1) . Simulations using more relevant 
progenitors would naturally support this comparison. 



6 CONCLUSIONS 

We have presented and discussed new visual wavelength-range 
observations of the innermost region of the supernova remnant 
SNR 0540-69.3 in the LMC. We have used the integral-field spec- 
trograph VIMOS/IFU that is mounted on the third unit telescope 
at the VLT of ESO in Chile. The IFU technique has up to now 
onl y been spar sely used with supernova remnants (cf., for exam- 
ple. lZharikov e t al. 2001). Our observations have provided us with 
three-dimensional maps of [Om]/15007 and [Sii] /1/16717, 6731 at 
an 0'.'33x0'.'33 spatial sampling and a velocity resolution of about 
31kms"'. We have decomposed the sulphur lines and created a 
"master" [S ii] line profile for each spectrum at high signal-to-noise 
ratio. We have also decomposed the oxygen lines to better probe 
high-velocity structures. In total we have about 1400 spectra with 
good [O III] and [S ii] line profiles. We have not only made intensity 
maps for [O m] and [S ii] , but also for the intensity ratio /[s „] //[o m] > 
and the average electron density. 

We have, moreover, recovered previous results of S05, M06, 
and LI 1 in terms of the "blob", line profiles, and overall projected 
structures on the sky. Additionally, with our VIMOS/IFU data we 
have obtained more detailed information. Assuming linear expan- 
sion of the ejecta, we have created a three-dimensional model, 
which reveals a structure that stretches from the position of the 
"blob" (Lll) and into the plane of the sky, at a position angle of 



PA^ 60°. Assuming a remnant age of 1000 years and an object dis- 
tance of 51 kpc, the structure has an inclination angle of about 65° 
relative to the line of sight. The position angle is close to the sym- 
metry axis of PA=45°-50°, which is sugested and discussed in LI 1 
to indicate present and past activities in the PWN. We speculate that 
the pulsar is positioned on this activity axis; in this case the pulsar 
has a velocity along the line of sight of a few hundred km s"' . 

The "blob" is most likely a region of shock activity, as it 
mainly emits in [Sii]. A decisive observation to test this specula- 
tion is to make velocity maps also of an ion of another element, viz. 
[O ii] /1/13726, 3729. The elongated structure of the central emis- 
sion is, furthermore, embedded in a volume of diffuse [O iii] emis- 
sion, which gives rise to a wide shoulder of a spatially integrated 
line profile for the full field. The shoulder is somewhat brighter 
towards bluer wavelengths compared to redder wavelengths. So, 
whilst ejecta mainly emit on the far side of the remnant, diffuse 
emission is biased towards the near side. 

We have, finally, compared SNR 0540-69.3 to other super- 
nova remnants and to a three-dimensional supernova model of 
(Hammer et al. ( 2010). In particular, the striking resemblance in X- 
rays between the PWN of SNR 0540-69.3 and the Crab - in com- 
bination with our findings in this paper - suggests that the sym- 
metr y axis is part of a torus in the PWN, which is in agreement 
with .Gotthelf & Wang ( 12000 ). and not a pulsar jet as was later sug- 
gested. Contrary to for the Crab, supernova ejecta are more affected 
by the torus activity in SNR 0540-69.3. The "blob" could be an un- 
usually active part of such an interaction. 
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